APCTP-Pre2008-002 
larXiv : 0806 . 4460i rhep-thl 
June, 2008 



00 
O 
O 



X 



Density Dependence of Transport Coefficients 
from Holographic Hydrodynamics 



Xian-Hui Ge*, Yoshinori Matsuo*, Fu-Wen Shu*, 
Sang- Jin Sin 1 "* and Takuya Tsukioka* 

* Asia Pacific Center for Theoretical Physics, 
Pohang, Gyeongbuk 790-784, Korea 

Qh' 

D \ gexh, ymatsuo, fwshu, tsukioka@apctp.org 

f Department of Physics, Hanyang University, Seoul 133-791, Korea 

sjsin@hanyang.ac.kr 

O 

\6 
o 

00 . 

, We study the transport coefficients of Quark-Gluon-Plasma in finite temperature 

and finite baryon density. We use AdS/QCD of charged AdS black hole background 
with bulk-filling branes identifying the £7(1) charge as the baryon number. We cal- 
culate the diffusion constant, the shear viscosity and the thermal conductivity to 
plot their density and temperature dependences. Hydrodynamic relations between 
those are shown to hold exactly. The diffusion constant and the shear viscosity are 
decreasing as a function of density for fixed total energy. For fixed temperature, the 
fluid becomes less diffusible and more viscous for larger baryon density. 



Abstract 



1 Introduction 



After the discovery of consistency on the ratio of the viscosity to the entropy density 
(i]/s) [1] in AdS/CFT correspondence [2-4] and RHIC (Relativistic Heavy Ion Collider) 
experiment, much attention has been drawn to the calculational scheme provided by string 
theory. Even some attempt has been made to map the entire process of RHIC experiment 
in terms of the gravity dual [5] . The way to include a chemical potential in the theory was 
figured out in the context of D4D8D8 setup [6,7]. Phases of these theories were discussed in 
D3/D7 setup and new phases were reported where instability due to the strong attraction 
is a feature [8-10]. 

Although QCD and M = 4 SYM are different, it is expected that some of the properties 
are shared by the two theories. It is an interesting question to ask how much one can learn 
by studying the various versions of AdS/CFT correspondence. The relevance is based on 
the universality of low energy physics. In this respect, the hydrodynamic limit is interesting 
since such limit can be shared by many theories in spite of the differences in UV limit. 

The calculation scheme of transport coefficients is to use Kubo formula, which gives 
a relation to the low energy limit of Wightman Green functions. In AdS/CFT corre- 
spondence, one calculate the retarded Green function which is related to the Wightman 
function by fluctuation-dissipation theorem. Such scheme has been developed in a series 
of papers [11-15]. 

For charged case, the calculations are more involved and corresponding works have 
been done only partially by various groups [16-19]. In [16,17], full thermodynamics for 
STU black hole [20,21] and the hydrodynamic calculations were performed for the (1, 0, 0) 
charge. In [18], charged AdS black hole background was considered and it was shown 
numerically that the ratio (rj/s) was l/(47r) with very good accuracy. Later, it was also 
proven that the ratio was universal in more general setup [19]. 

In this paper, we perform hydrodynamic calculation directly for Reissner-Nordstrom- 
Anti-deSitter (RN-AdS) black hole, which corresponds to the (1, 1, 1) charged STU black 
hole. Master equations for the decoupled modes are worked out explicitly. Vector modes 
of Maxwell field (as well as the vector modes of metric) have diffusion pole, contrary to 
the chargeless case. As a consequence, the diffusion constant is calculated and the ratio 
(rj/s) = l/(47r) is shown to be exact in this case also. We observe that the diffusion 
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constant and the shear viscosity decrease as we increase the charge with fixed total system 
energy (or equivalently the black hole mass), while the shear viscosity increases for fixed 
temperature. 

The charge in RN-AdS black hole is usually regarded as -R-charge of SUSY [25]. We 
here consider an another interpretation in the following way: One can introduce quarks and 
mesons by considering the bulk-filling branes in AdSs space. The overall U(l) of the flavor 
branes is identified as the baryon charge. The U(l) charge in this model [26] minimally 
couples to the bulk gravity since the bulk and the world volume of brane are identified. 
RN-AdS metric can be considered as the consequence of the back reaction of the AdS 
black hole to the charge. Therefore the U(l) charge in RN-AdS can be identical to the 
baryon charge. As a result, we can calculate the transport coefficients in the presence of 
the baryon density. 

One can give an explanation of hydrodynamic mode in meson physics. In our interpre- 
tation, the Maxwell fields are the fluctuations of bulk-filling branes, therefore they should 
be interpreted as master fields of the mesons. Then hydrodynamic modes are lowest lying 
massless meson spectrum. In terms of brane embedding picture, this massless-ness is due 
to the touching of the brane on the black hole horizon. Near the horizon, the tension of 
the brane is zero due to the metric factor and it can lead to the massless fluctuation. Then 
the massless spectrum can not go far from the horizon in radial direction. In this picture, 
hydrodynamic nature is closely related to the near horizon behavior of the branes. 

This paper is organized as follows: In section [21 we introduce RN-AdS black hole and 
review correlation function calculation at finite temperature in AdS/CFT correspondence. 
In section [3], a formulation on the metric and the gauge perturbations in RN-AdS back- 
ground is given. We then calculate the correlators in hydrodynamic regime and obtain 
the diffusion pole in section HI In section [5J the shear viscosity is calculated via Kubo 
formula. We also show that the result is consistent with the hydrodynamic relation of 
diffusion constant and the viscosity. The thermal conductivity is also calculated in this 
section. Conclusions and discussions are given in section EJ Three appendices are given to 
provide the details of the calculations. 
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2 Basic Setup 

2.1 Minkowskian correlators in AdS/CFT correspondence 

Before introducing RN-AdS black hole, we briefly summarize Minkowskian correlators in 
AdS/CFT correspondence. We follow the prescription proposed in [11]. Let us consider 
fluctuations of fields which satisfy equations of motion at the linearized order. We work 
on the five-dimensional background, 

ds 2 = g^dafdx" + g uu {du) 2 , (2.1) 

where and u are the four-dimensional and the radial coordinates, respectively. We refer 
the boundary as u = and the horizon as u — 1. A solution of the equation of motion 
may be given, 

<P(u,x)=J-0^e^f k (u)Mk), (2.2) 

where fk(u) is normalized such that fk(0) = 1 at the boundary. An on-shell action might 
be reduced to surface terms by using the equation of motion, 

s[<M = J ^^(-*O0(M)0o(*O[^- (2-3) 

Here, the function Q(k,u) can be written in terms of f±k{u) and d u f±k(u), for example, 
for a scalar field, 

G(k,u) = K^g-g uu U(u)d u f k (u), 

with some constant K. The direct generalization of AdS/CFT correspondence, or Gubser- 
Klebanov-Polyakov/Witten relation [3,4], to Minkowski spacetime gives the relation, 

e iJ> eA = e iS[4>o] 



where the operator O is defined in the boundary field theory. From this relation, one 
may obtain a Green function by taking second derivative of the action with respect to the 
boundary value of the field, 



G(k) = —G(k,u) -Q{-k,u) 

u=0 



u=l 



u=l 



u=0 

However, this quantity is real and cannot be a retarded Green function. This can be seen 
as follows. The imaginary part of Q(k,u) is proportional to a conserved flux. Then, its 



contributions at the boundary u = and at the horizon u — 1 cancel completely. Even if 
one neglects the contribution from the horizon, G(k) is still real. The reality condition of the 
equation of motion implies Q(—k,u) = Q*(k,u), and the imaginary part of G(k) vanishes 
again. Therefore we should impose the "retarded" condition to the Green function. 
Son and Starinets proposed that the retarded (advanced) Green function is given by 



with incoming (outgoing) boundary condition at the horizon. Generally, the contribution 
at the horizon is oscillating and averaged out to zero. In order to avoid this, we have 
to consider incoming or outgoing boundary condition. Taking away the contribution at 
the horizon, we obtain Q(k,u) with a non-zero imaginary part. Physics at the horizon 
affects the Green function only through the boundary condition. In general, there are 
several fields in the model. We write the Green function as Gij(k), where indices i and j 
distinguish these fields. The surface terms are always associated with equations of motion. 
We choose the former index to indicate the field whose equation of motion is associated 
with the Green function. 

In this paper, we work in RN-AdS background and consider its perturbations so that 
essential ingredients are perturbed metric field and U(l) gauge field. Here we define the 
precise form of the retarded Green function which we discuss later: 



where the operators T lxu (x) and J M (x) are energy-momentum tensor and U(l) current which 
couple to the metric and the gauge fields, respectively. 

2.2 Reissner-Nordstrom-AdS background 

In this paper, we consider N c D3-branes and Nj D7-branes, and treat the D3-branes as a 
gravitational background. The D7-branes are wrapping on S* 3 of S 5 , and we neglect this 
S 3 dependence. We do not consider the perpendicular fluctuations of D7-branes, and the 
effective action then becomes that for five-dimensional gauge theory. If the D7-branes touch 



G R 



(k) = 2g(k,u) 



(2.4) 



u=0 




(2.5) 
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l - J d 5 xV=g- Tr {T mn T mn ), (2.6) 



the D3-branes, the D7-branes fill the AdSs completely. The induced metric on the D7- 
brane is identical to the AdS bulk metric [26] . This model corresponds to A^ = 4 SYM with 
massless quarks. If we introduce the baryon charge at the boundary theory, its chemical 
potential is identified as the tail of the U{\) gauge potential on the flavor brane [6,7]. (See 
also [8-10,27,28] for later development.) We consider the phenomenological model taking 
only AdSs part and neglecting S 5 part. Then there is no way to distinguish the bulk gauge 
field and the brane field. The baryon charge and the .R-charge have the same description 
in terms of the U(l) gauge field living in the AdSs space. A charged black hole (RN-AdS 
black hole) is then induced by its back reaction. This corresponds to the M = 4 SYM in 
finite temperature with finite baryon density. It is the case that we consider in this paper. 

The action for the gauge field dual to the baryon current is given by the U(l) part of 
the Dirac-Born-Infeld actiorj* 

SW - 



4e 2 

where the gauge coupling constant e is given by 

I _ N c Nf 

J ~ (2^' ( ' 

with I the radius of the AdS space. Notice that the gauge field is that of the diagonal U(l) 
of U(Nf) flavor brane dynamics, which is dual to the baryon current at the boundary. 
Together with the gravitation part, we arrive at the following action which is our starting 
point: 

S[g mn , A m] = ^J d 5 x (r-2K)--^J $x^-gT mn T mn , (2.8) 

where we denote the gravitation constant and the cosmological constant as k 2 = 87TG5 and 

A, respectively. The U(l) gauge field strength is given by T mn {x) = d m A n (x) — d n A m (x). 

The gravitation constant is related to the gauge theory quantities by 

I 3 N 2 r , . 

^ = 4^- 

Suppose we have baryon charge Q. This should be identified to the source of U{\) charge 
on the brane hence on the bulk. Then we can relate it to the parameter in RN black hole 
solution by considering the full solution to the equation of motion, 

Rmn ~^9mnR fi'mnA ^ -^mn> (2.10) 



* The indices m and n run through five-dimensional spacetime while /x and v would be reserved for 
four-dimensional Minkowski spacetime. Their spatial coordinates are labeled by i and j. 



where energy-momentum tensor T mn (x) is given by 

An equation of motion for the gauge field A m (x) gives Maxwell equation, 

V m F mn = ^L^fvqfV't^i, - diAfSj = 0. (2.12) 

Here we assumed that there is no electromagnetic source outside the black hole. One can 
confirm that the following metric and gauge potential satisfy the equations of motion (I2.10p 
and 

d* 2 = £( - /(r)(dt) 2 + X>*) 2 ) + ^4t(^) 2 , (2.13a) 



i=l 



A t = -%+^ (2.13b) 



with 



ml 2 q 2 l 2 6 
if and only if q is related to the Q by 



= (2,4) 

It should be noted that a ratio of the gauge coupling constant e 2 to the gravitation constant 
k 2 is 

e 2 N c ~ . 

- = w l 2 - 15 

K z JMf 

Since the gauge potential At(x) must vanish at the horizon, the charge Q and the chemical 
potential fi are related. The parameters m and q are the mass and charge of AdS space, 
respectively. This is nothing but Reissner-Nordstrom-Anti-deSitter (RN-AdS) background 
in which we are interested throughout this paper. 

The horizons of RN-AdS black hole are located at the zero for /(r)Q, 

lir) = 1 - ^ + ^ = ^ (r 2 - rl) (r 2 - ri) (j* - $ , (2.16) 



In order to define the horizon, the charge q must satisfy a relation g 4 < 4m 3 / 2 /27. 
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where their explicit forms of the horizon radiuses are given by 

^ = ($( 1+2cos (5 + r))) (2 - 17a) 



/•_ I I I ^2 cos 



r n = 



' (J)))", ^ P.17B, 

(i( 1+2Ms (i + f))) '• < 2 - 17c > 



with 



= arctan — — — — ; 

V 2m 3 / 2 - 27g 4 / 



and satisfy a relation + r 2 = — r\. The positions expressed by r + and r_ correspond to 
the outer and the inner horizon, respectively. It will be useful to notice that the charge q 
can be expressed in terms of 9 and m by 



4 4m 3 / 2 . 2 /0 



o = sin 

H 27 V 2 



The outer horizon takes a value in 



/ in 

^j-l <r 2 + < V^l, 



where the upper bound and the lower bound correspond to the case for q = and the 
extremal case, respectively. 

We shall give various thermodynamic quantities of RN-AdS black hole [25,26]. The 
temperature is defined from the conical singularity free condition around the horizon r + , 

r%f'(r+) t +(a \q 2 ^ 

+ 

where we defined the parameters a and b by 

q 2 l 2 , I 2 
r% 2r + 

In the limit q — > 0, these parameters go to 

/3/2 

and the temperature becomes 



a = V> & = — • ( 2 - 19 ) 



m. 



1/4 



The entropy density s, the energy density e, the pressure p, the chemical potential fi and 
the density of physical charge p can be also computed as 

(2.21) 

1 + aV (2.22) 



s = 


2ixr\ 


Til 3 




An 2 b 3 ' 


e = 


3m 


3^ 3 


2k 2 / 3 
e 


32k 2 6 4 


P = 


3' 




fi = 


Q 

r\ ' 




P = 


2Q 
e 2 / 3 ' 





(2.23) 
(2.24) 

(2.25) 



3 Perturbations in RN-AdS Background 



In RN-AdS background, we study small perturbations of the metric g mn (x) and the gauge 
field A m (x), 

a = o (0) + h 



A = -i- A 



(3.1) 
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where the background metric gm n (x) and the background gauge field A™ (i) are given 
in (I2.13al) and (|2.13bl) . respectively. In the metric perturbation, one can define a inverse 



metric as 

g mn = g^ mn - h mn + 0(h 2 ), 

and raise and lower indices by using the background metric gmn(%) and g^ mn (x). 

Now we shall consider a linearized theory of the symmetric tensor field h mn (x) and the 
vector field A m {x) propagating in RN-AdS background. In the first order of h mn (x) and 
A m (x), the Einstein equation (12.101) can be written as 

_ !„(0) d(1) _ \ d(0) , u A _ 2 T (1) /oo\ 
ll mn 2^ mn 2 mn IL mn 1 ^ — rv J- mn - \ -^J 

In the expression above, the scalar curvature (x) is constructed by using the background 
metric gm n {x) and the following tensors are newly defined: 



Rmn ~ 2 (^kV m h n k + VkV n h m k — VkV k h mn — V m V T 



h 
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T (l) _ 1 ( _ p(0) p(0) h kl , 1(0) p(0) p(0) p h kl _ \ F (0) ^(0) 
J mtl _ e 2 \ r mk r nl 11 ' 2^ mn k P 1 ^ hmnF kl r 

' r mk r n ' 1 nk r m n 9mn £ kl r 



mk*- n ^ ± nk ± c^ran^ kl 

where the Ricci tensor Rmn(x), the covariant derivative and the field strength Fm n (x) 
are defined through the background metric gmn{x) and the gauge field Am(x). We de- 
note a trace part of the metric and a field strength for the perturbative parts as h(x) = 
h mn g^ mn (x) and F mn (x) = d m A n (x) — d n A m (x), respectively. On the other hand, the 
Maxwell equation (I2.12p becomes 



= V m (F m ™ - F {0)m k h nk + F {0)n k h mk + ^F (0)mn h^ 
= - ? i=a m {v / ^^ (0)m ^ (0)ni (9 fc ^ - d t A k ) 

_ F (0)m kh nk + F (0)n kh mk + _ F (0)mn^ j _ (3 3) 

The above equations of motion (I3.2p and (13. 3p can be derived from the following action: 
S\h m ,A m ] =~ J d s x^mh m h"" l VJ, - V m h« k V„h m k 



+ 



^L F (0) F (°)hmk h 
g2 r mn r H 11 " 



2 



nl 



-2{F^ k F n k h mn + F^Fjh™ - ±F®F™h)y 

(3.4) 

By using the equations of motion, an on-shell action is reduced to surface term 
S[ht n , At] = ~ J d 5 a;9 m | (h mn V n h + h\/ n h mn - 2h nk V n h m k 



+h kl V m h kl - hV m h 



.±Jd 5 xdJV^gV)A n (F« 
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F (0)m kh nk + F (0)n kh mk + _p{0)van h 



(3.5) 



We shall work in the h rm {x) = and A r (x) = gauges and use the Fourier decompo- 



where we choose the momenta which are along the z-direction. In this case, one can 
categorize the metric perturbations to the following three types by using the spin under 
the 0(2) rotation in (x, y) -plane [12]: 

• vector type: h xt ^ 0, h xz ^ 0, (others) = 

^equivalently, h yt ^ 0, h yz ^ 0, (others) = 

• tensor type: h xy ^ 0, h xx = —h yy ^ 0, (others) = 

• scalar type: h tz ^ 0, h tt ^ 0, h xx = h yy ^ 0, and h zz ^ 0, (others) = 

We consider the first two types in this paper. The scalar type perturbation would be 
studied elsewhere. 

3.1 Vector type perturbation 

In this subsection, we study the vector type perturbation in RN-AdS background. From 
explicit calculation, one can show that only x-component of the gauge field A x (x) could 
participate in the linealized perturbative equations of motion. Thus independent variables 



h xt {x) + 0, h xz {x) ^ 0, A x {x) ^ 0, (others) = 0. 

We start by introducing new field variables, h x (r) = g^ xx h xt (r) = (I 2 / r 2 )h xt {r) and 
h x {r) = g^ xx h xz (r) = (I 2 /r 2 )h xz {r). Nontrivial equations in the Einstein equation (13. 2ft 
appear from (t,x), (r,x) and (x,z) components, respectively: 



sition 




are 




(3.6a) 
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Qq 2 l 2 u 



= kfhf + uhf + -^^A x , (3.6b) 



Qr 

( r 5f)> ;4 

g I rtl \ J ) LI' , 



K" + ^lfK> + — (ukhf + u?K) , (3.6c) 



r 5j z r 4 / 5 

where the prime implies the derivative with respect to r. In the set of equations, the 
equations (I3.6al) and (13.6bj) imply (I3.6cp . On the other hand, in the Maxwell equation 
(13. 3p . the x-component gives a nontrivial contribution, 

° = k + { 4ff A '* + ^ 2 - k2 fh + (x7) 

Taking the limit in which the charge q goes to zero, the metric and the gauge perturbations 
are completely decoupled. 

We now look for solutions of our set of equations. First of all, from the equations (I3.6ap 
and (I3.6bp . we can obtain a second order differential equation for hf{r) and A x {r), 

= hf" + <£QV + ^ ( B (rV)' + ^ 2 - k 2 f) ) hi 
, 6g 2 / 2 /^ , (r- 1 /)' A , , Z 4 " 2 



Q \r 5 r 4 / x r 9 / 2 
Together with the equation of motion (13. 7p . we treat hf(r) and Aj.(r) as independent 
variables. Having the solutions for these, one can get one for h x J{r) by using the equation 
(13.6bp . In order to solve these equations, we find it is useful to introduce linear combinations 
of the variables 

80 4 5ux/ ( 3a/ 4 C± 2 \ . 
$± = — 7T r h t + ( ~ 77TT7 + "7T r ) A x , (3-9) 



with constants 



C± = (1 + a) ± + a) 2 + 3ab 2 k 2 , 



so that we can obtain second order ordinary differential equations in terms of these new 
variables. In fact, the equations of motion ( 13. TP and ( 13.8ft could be rearranged as 

= * ± " + *±' + ^i(- 2 - * 2 /)* ± - |§j* ± . (3.10) 

In the chargeless limit, the two equations of motion (13.101) for $+(r) and $_(r) give de- 
coupled ones for A x (r) and hf(r), respectively. 

We will consider these equations of motion in low frequency limit so-called hydrody- 
namic regime. In the hydrodynamic regime we could obtain the diffusion pole and the 
thermal conductivity from retarded Green functions. 
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3.2 Tensor type perturbation 

Next we shall focus on the tensor type perturbation. By considering the spin or by cal- 
culating directly, the metric perturbation is decoupled from the gauge perturbation. Thus 
independent variables are 

h xy (x) ^ 0, h xx (x) = -hyy(x), (others) = 0. 

A nontrivial equation of motion in ( 13. 2p is coming from (x,y) component. As we did in 
the vector type perturbation, it might be convenient to introduce new variable hy(r) = 
g^ xx h xy (r) = {I 2 / r 2 )h xy {r) . We then get the following equation of motion: 



An another equation of motion for h xx (r) = —h yy (r) is as the same form of ( 13.111) . We use 
this equation of motion to study the shear viscosity in the hydrodynamic approximation. 

The equations ( 13.101) and ( 13.111) can be rewritten as Schrodinger-like equations through 
suitable field redefinitions. Their potentials were derived by Kodama and Ishibashi [29]. 

4 Diffusion Pole in Hydrodynamic Regime 

In the hydrodynamic regime, it is standard to introduce new dimensionless coordinate 

u = r\/r 2 which is normalized by the outer horizon. In this coordinate system, the 

horizon and the boundary are located at u = 1 and u = 0, respectively. Defining the new 
A (u) l A 

variable B(u) = — = -A x (u) where a is the chemical potential given by (12.241) . 

u 4Qo^ 

our basic equations (I3.6al) - ()3.6cp and (13. 7p are rewritten in this new coordinate system: 




(3.11) 



I b 2 / 
= hf - -hf (ukh* + k 2 h 



u uf\ 
= kfhf + uhf - 3auuB, 



X 



- 3auB 



(4.1a) 



(4.1b) 




(4.1d) 



(4.1c) 



with 



f(u) = (1 -u)(l + u-au 2 ). 
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Here the prime now means the derivative with respect to u. The equation (13.101) may be 
also written down as 

o = $±" + ^ff-**! + k2 f) $ ± - % $ ±> ( 4 - 2a ) 

u 2 f uf 2 \ J f 

for 

$± = -hf - 3aB + —B. (4.2b) 

u u 

Getting the solution for $ ± (ii), one can access to solutions for hf(u) and B(u), 

hf = «$_ + ^V« 2 (*+-*-)- ^T^ M ( $+ " $ -) ' (4 ' 3a) 

5 = 1 M C$ + - $_) . (4.3b) 



The constants C± could be expanded in this regime, 

C + = 2(1 + a) + ^^yfc 2 + 0(A; 4 ) 
3n/) 2 

2(1 + a) V 7 



(4.4) 



First, let us consider the equation for $_(it). Following the usual way to solve differ- 
ential equations, we impose a solution as $_(«) = (1 — u) u F^(u) where F_{u) is a regular 
function at the horizon u = 1. Substituting this form into the equation of motion, one can 
fix the parameter v as v = ±iuj / (4-7rT) where T is the temperature defined by the equation 
( 12.181) . We here choose 



v = —i- 



'4ttT' 

as the incoming wave condition. 

Now we are in the position to solve the equation of motion in the hydrodynamic regime. 
We start by introducing the following series expansion with respect to small uj and k: 

F_(u) = F (u) + wFi(u) + ^diu) + 0(u 2 , uk 2 ), (4.5) 

where F (u), Fi(u) and G\(u) are determined by imposing suitable boundary conditions. 
In order to do the perturbative analysis, it might be convenient to rewrite the equation 
fOajl for $_(«) as, 



= (u 2 (l - u){l +u- au 2 )FJ^j 
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+iuj^—u 2 (l + u- au 2 )F'_ + iu—^—u(2 + 3u - 4aw 2 )i+ 
2 — a 2 — a 



u 



(2 - a) 2 1 + u - au 2 
x ((2 - a) 2 + (1 - a) (3 - a)w + (1 - 4a + a 2 )u 2 - a(2 - a)w 3 + aV)i+ 



-ifc 2 &Vl 



3 a 



2(1 + a) 

The solution can be then obtained recursively. The result is as follows: 



(4.6) 



Fq(u) = C, (const.), 



/•■(„) = }Ch{ - \ + 2a / ^ —{ leg 



2 VI + 4a(2 - a) 



+ M + 2(2-a) 
Cb 2 f 1 
2(1 + a) V + m 



Vi - 


1- 4a 


1 - 


2a 



/ l-2au\ 
/ l + . \ 



(4.7a) 



- log 



log 



V 1 vT+~4a/ 

l + u — au 2 
2 - a 



Vi- 


f- 4a 


l - 


2a 



(4.7b) 
(4.7c) 



All of the solutions should be regular at the horizon u — 1 and the functions Fi(u) and 
G\(u) should be vanished there. The constant of integration C will be estimated later. 

Next, we shall study the equation for $+(«)• It might be useful to introduce new 
variable 



3a 



2(1 + a) U/ 

In terms of new variable, the equation of motion (j4.2aj) for becomes 



(4.8) 



o = $" + 



((- 



3a \2 x' 



2(1 + a 



3a \ 2 
u) f 



2(1 + a 



u/ 5 



^ 2 -fc 2 / 1 



3 a 



2(1 + a) 



u $ 



(4.9) 



Assuming again <&+(«) = (1 — u) u F{u) where is a regular function at u — 1, the 

singularity might be extracted. The equation of motion (14. 9p becomes 







1 — u I ( 1 + u — au' 



3a u) 2 f" 



2(1 + a 



+2iu 



2- a 



1 + u — au 111 — 



3a \ 2 ~, 
u ) F 



2(1 + a 



The derivation of the solutions is given in Appendix A. 
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+IUJ 



2- a 



2h2 



1 + u — au 111 



3a uYXf 



2(1 + a) 



u 2 b 



__ 3a a 2 



(2 -a) 2 w(l + w-aw 2 ) 

x ( (2 - a) 2 + (1 - a) (3 - a)u + (1 - 4a + a 2 )u 2 - a(2 - a)u 3 + au 4 )F 



k 2 b 2 ( 3a w 3a \ 2 ~ 
1 + ZT7Z ru) 1-— r« F, 



it \ 2(1 + a) A 2(1 + a) / 

where we used the incoming wave condition v = —iu/(4nT) as same as before. 
We impose a perturbative solution as 

F{u) = F Q {u) + wFi(u) + fc 2 C?i(u) + 0(u 2 , uk 2 ), 

and then we obtain the following result^: 

Fq(u) = C, (const.), 

(76 f 27a 2 



1 -w 



2-all + 4a\2 + 2a- 3aw 



1 - 10a - 2a 2 . 
+ — . . ^ ( log 



G x (u) = CJ[u 
= Cb 2 



2(1 + 4a) 3 / 2 



1 , / 1 + u — au 2 
H — log 

2 & V 2-a 



9a 2 (14 + 31a + 8a 2 



/ 1 - 2au \ 

VI + 4a 
1 -2a 



log 



/ l-2au\ 
1 + ; X 
VI +4a 

1 - 2a 



V yrrw 



1 - it 



2(1 + a)(l + 4a) (2 - a) 2 \2 + 2a - 3aw 
(1 + a) ^3a(2 - a)(5 + 2a) - 2(1 + a)(l - 10a - 2a 2 ) log(3a; 

(2-a) 3 (l + 4a) 3 / 2 



x lot 



/ 1 - 2aw \ 

VI +4a 
1 -2a 

V ~ yrrs/ 



- log 



I 1 + . \ 

yi + 4a 



1 - 2a 



. 1+ . 

V yrT4^/ 



4(1 + a)- 
' (2 - a) 3 



log u log (l — u 



(4.10) 



(4.11) 



(4.12a) 



(4.12b) 



' The detail is given in Appendix B. 
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(1 + a) (9a{2 - a) + 2(1 + a)(l + 4a) log(3a)) 

(2-a) 3 (l + 4a) 
54a 2 (l + a) ^ wlogu 



los 



1 + it — aw z 



(2 - a) 2 (l + 4a) \2 + 2a - 3au 
4(1 + a) 



7T 



L12H - y 



(2 - a)3 

2(1 + a) 2 
(l + 4a) 3 / 2 (2-a) 3 

x ^1 - 10a - 2a 2 + (1 + 
x ( log u log I 1 



2au 



1 - Vl + 4a 



A- 



+Li 2 



2aw 



1 - VI + 4a 
(l - 10a - 2a 2 - (1 +4a) 3/2 



Li 2 



+ log(3a) log 



2a 



2au 



1 - VTT4a" 



2a 



\ 1 - v/TTIa/ 



1 - VI + 4a 



x ( log u log ( 1 



2au 



1 + VI + 4a 



+Li 2 



2au 



Lig 



+ log(3a) log 



2a 

1 + VI + 4a 



2au 



\ 



1 + VTT4^ 



2a 



\ 1 + vTTW 

(4.12c) 



1 + VI + 4a, 

where Li 2 (u) is the polylogalithrrlj]. It should be mentioned that the defined functions H(u) 
and J(u) in Fi(u) and G\{u) are finite at the boundary u — 0. 

Let us consider the integration constants C and C . These could be estimated in terms 
of boundary values of the fields 

lim^H = (^)°, limh:(u) = (h:)°, limS(«) = 

u^O it— >U M— >U 



Taking a derivative of $±(it) and using the equation of motion (14. lap , we can get relations 

u 2 & ± - C±uB' = — Uukh x + k 2 h x ) - C±B. 



* The polylogalithm appears from 



(Li 3 («))' 



log(l - it) 



Some values are given as, Lia(— 1) = — 7r 2 /12, Li2(0) = and Li2(l) = 7r 2 /6. 
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We evaluate the equations above at the boundary, 

lim(«V ± - C ± uB'\ = b 2 (uok(h x )° + k 2 {h x t f\ - C±(B)°, (4.13) 

so that we may fix the constants C and C from + parts, respectively, 

3ab 



C = r , (4.14a) 

7 2 

2(T+a) 

-b 2 U{Kf + k\Ktf) + (2(1 + a) + ^rk 2 ) (B)° 

C = = [ +a) , (4.14b) 

l + cuH(0) + k 2 J(0) 

where we used the obtained solutions $±(m) and the relation (14.3bl) for B'{u). It should 
be noted that the boundary value of uB'{u) is vanished. In the equation ( I4.14al) . one can 
see the existence of the hydrodynamic pole in the complex cu-plane. 

Now we proceed to calculate the Minkowskian correlators. For the vector type pertur- 
bation, the on-shell action (13.51) becomes 

S[h' t ,K,B\ = f ^Ll-h' t (-k,u)Kfik,u) - \ht(-k,u)ht(k,u) 

• -h x (-k,u)h x '(k,u) + i ^Y-h x (-k,u)h x (k,u) 



u u 

u=l 



-3af(u)B(-k, u) (B'(k, u) - J^ h t( k , « 



u=0 

(4.15) 



Using the obtained solutions, we can lead the following relations between the radial deriva- 
tive of the fields and their boundary values near the boundary u = e: 

h x '(e) = -b 2 Lk(h x z )° + k 2 (h x ) 



e 

+ 



b 

iuj 



-k 2 



b(uk(h x z f + k 2 (h x f^ + 3iauj(B)° + 0(u 2 k, ujk 2 )\ 



2(1 + a) 

+0(e 2 ), (4.16a) 
h x z '(e) = b 2 {u; 2 (h x z ) + u J k(h x f 

S \ - f n / , ™. \ n i / 1 fv \ n \ 3 ft & 



b(cu 2 (h x )° + iuk(h x ) ) + uk(B)° + 0(uj 2 k, uk 2 

V / 2(1 + a) 



2(1 + a) 
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B'(e) 



+0(e 2 ), (4.16b) 



& ,2 1 2 



; rk 

2(1 + a) 



+ (o 2 £; 2 (5) + 0{uk 2 )^) loge + 0(e). (4.16c) 

By using the relation (12.41) and the definition (12.51) . we can read off the correlators in the 
hydrodynamic approximation, 

I 3 ( k 2 



G — ( ^ ) = iwU^»J' (4 - 17a) 

G xt X M k) = Gxz k) = - (-^) , (4.17b) 

I 3 ( LU 2 



G xt Jp, k) = G x xt (co, k) = ( - lU — ) , (4.17d) 

e A l 6 \uo — Dk l J 

_ , 3a/ / iu \ (2 -afl 

° x ^ k) = 4(l + a)frV ^TZT^F J " 8(iTW lW ' ( } 



where we subtracted the contact terms. In the final expression above we rescaled the gauge 

AQb 2 

field (B)° to the original one (A x )° = ^ (B)° and raised and lowered the indices by using 
the flat Minkowski metric rj^ = diag(— , +,+,+) in the four-dimensional boundary theory. 
Taking the limit in which the charge q goes to zero, the results coincide with the known 
ones in [12]. In this limit, the correlators (14.17dj) and (I4.17el) vanish, while the correlator 



( I4.17fj) has no diffusion pole and the subleading term reproduces the consistent result. The 
same interesting structure was found in the single (1,0,0) -R-charged black hole [17]. The 
constant D is the diffusion constant 

with _ 

9 ~ arctan ( ^W^l 2 -27 q * \ 
V 2m 3 / 2 - 27g 4 J 
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All of the correlators in the vector type perturbation exhibit a diffusion pole. The behavior 
of the diffusion constant is drawn as a function of the charge q and the mass m in Figure 
[H and as a function of the charge q and the temperature T in Figure EJ In the chargeless 




Figure 1: D vs. q and m (I = 1) Figure 2: D vs. q and T (I = 1) 

limit, the diffusion constant becomes 

D^D 

where the temperature T is given in (12.201) . 

5 Shear Viscosity in Hydro dynamic Regime 

In this section, we solve the equation of motion (13.111) in the hydrodynamic regime and 
obtain the shear viscosity. We could also see the hydrodynamic relation and the formulation 
of the thermal conductivity. 

After changing the coordinate r to u = r^/r 2 , the equation (13.111) can be rewritten as 

with 

f(u) = (1 -u)(l + u-au 2 ), 



1 
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where the prime means the derivative with respect to u. Removing the singularity around 
u — 1, the equation becomes 



= (-(1 - + u- au 2 )F'^j 



-\ , 2l> -(1 -«)(! + «- au 2 )F' - «j -^-(1 + au 2 )F 



(2 - a) u 



(2 - a) ti 2 



(2-a)V(l + w-aw 2 ) 

x ( (a - 2) 2 + (a - 3)(a - l)u + (a 2 - 4a + l)w 2 + a(a - 2)w 3 + aV ) F 



-A: 2 -F, 
tr 

where we imposed the incoming wave condition 

h x y ( u ) = (i- u y iuj/ ^ nT) F( u ). 

Perturbative solutions for F(u), 

F{u) = F (u) + uFi(u) + A; 2 Gi(u) + O(oo 2 , ook 2 

can be obtained aj*| 

Fq{u) = C, (const.), 
F x {u) = CH{u) 

3 



(5.2) 



(5.3) 



(5.4) 



(5.5a) 



Cb 



2(2 - a) \ 

+ log 



log 



/ 1 - 2au \ 
y/1 + 4a 



1 



1 - 2a 



-log 



/ 1 - 2au \ 
Vl + 4a 



1 + u — au 2 



1 - 2a 



V 1 ' VTTW 



2 - a 



Gi(u) = CJ(u) 
Cb 2 



log 



/ l-2aw\ 
VI + 4a 



1 



1 - 2a 



log 



/ 1 - 2au \ 

I 1+ ; \ 

VI + 4a 



1 - 2a 



. 1+ , , 



(5.5b) 



(5.5c) 



\ VTT4a7 

Since the function hy(u) goes to (h*) at the boundary u — 0, the constant C can be fixed 

as 



C 



x\0 



l + uH(0) + k 2 J(0)' 



(5.6) 



The detail is given in Appendix C. 
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Taking the limit q — > 0, the solution recovers the result in [12]. The solution of h x (u) is 
the same form as hy(u). 

Let us evaluate the Minkowskian correlators. The relevant part of the metric pertur- 
bation in the on-shell action (13. 5p becomes 



+^h x J-k,u)h x x '(k,u) - ^h x J-k,u)h x x (k,u) 
u u z 



u=l 

(5.7) 



Near the boundary u = e, using the perturbative solution for h x (u), we can obtain 

h x/ (e) = eb(ico + bk 2 ^ (h x )° - b 2 k 2 (h x )° + 0(u 2 , uk 2 ). (5.8) 

The same relation for h x (u) might be satisfied. Therefore we can read off the correlation 
functions from the on-shell action (15.71) . 

G X y xyiyi i k) = G xx xxi^i k) = Gyy yy(uj, k) 
73 , v 

iu + bk 2 ), (5.9) 



16k 2 6 3 

where we subtract contact terms. 

The result above can be used to estimate the shear viscosity 77 via Kubo formula, 

77 = - km = . (5.10) 

00 loKrtr 

Therefore we can conclude the following relation between the shear viscosity 77 and the 
entropy density s which is given in the equation (12. 21ft : 

2=1. (5.11) 
s Air 

The behavior of the shear viscosity is drawn as a function of the charge q and the mass m 
in Figure [3] and as a function of the charge q and the temperature T in Figure HI 

In hydrodynamics, the following relation is held: 

D = ^- } (5.12) 
e + p 

where e and p are the energy density and the pressure defined in ( 12. 22ft and (12. 23f) , respec- 
tively. Using the obtained diffusion constant (I4.18p . the shear viscosity could be calculated. 
We can confirm the result coincides with fl 5 . 1 [) which was obtained from Kubo formula. 
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Figure 3: r] vs. q and m (k = I = 1) 



Figure 4: 77 vs. q and T (k = I = 1) 




The thermal conductivity can be also computed from the Green function by using 
Kubo formula [17], 

^JL^ MM , (5 , 13) 

where the density of physical charge p is given by (12.251) . Here we can use the retarded 
Green function G x x (uj, 0) given by (14.17fj) as G(u, 0). Thus we obtain 

k t = 2vr 2 — -L. = 27T 2 -^^. 5.14 
y K J /i iV/ p 2 

The behavior of the thermal conductivity kt is drawn as a function of the charge q and 
the mass in Figure [5] and as a function of the charge q and the temperature T in Figure [H 

6 Conclusions and Discussions 

In this paper we considered holographic QCD in the presence of the baryon density by 
introducing the bulk-filling branes. We use RN-AdS black hole geometry as the gravity 
dual of such system. We have seen the diffusion pole structure in vector type perturbation. 
It is worth mentioning that the correlator of Maxwell fields in the vector mode G x x (u, k) 
has the diffusion pole unlike the charge free case. The transport coefficients have been 
calculated in holographic hydrodynamics and their temperature and density dependence 
was demonstrated. 

The diffusion constant decreases as charge increases for fixed temperature. Physically, 
this implies that the fluid is less diffusible for large baryon density. By calculating the shear 
viscosity analytically, we showed that the shear viscosity rj and the entropy density s satisfy 
the universal ratio (r]/s) = l/(47r) which has been originally suggested in [1]. For fixed 
temperature, the fluid becomes thicker as the charge increases. We have also seen that the 
diffusion constant and the shear viscosity satisfy the suitable relation for hydrodynamics. 
The calculation of the thermal conductivity shows that it satisfies (an analogue of) the 
Wiedemann-Franz low. 

It is very interesting to study the pole structure of scalar type as well as vector type 
of gravitational perturbations. Also it is important to carry out higher order calculations. 
Such result will be useful to get the higher order transport coefficients in the presence of 
the conserved current. We will report on these issues in the forthcoming publications. 

In our interpretation, the fluctuations of bulk-filling branes are regarded as master fields 
of the mesons. Near the horizon, the tension of the brane is zero due to the metric factor 
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and it can lead to the long range fluctuation. This becomes the hydrodynamic mode. One 
important question is about the meaning of hydrodynamic mode in terms of meson physics. 

Further question in this direction is how we can understand the dispersion relations of 
vector modes of Maxwell fields in terms of the particle spectrum with dissipation. This 
vector mode cannot propagate in neutral medium while it can in charged medium. In 
addition, the tensor mode does not propagate in the medium. It is interesting to consider 
their interpretations in terms of meson physics. More thought on these points is to be 
pursued in the future. 
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Appendix A. Perturbative solutions for <£>_ 

From the equation of motion (14. 6p . one can read off one for F (u), 



= (u 2 (l - w)(l +u- au 2 )F^j . 



(A.l) 



A general solution is given by 
F (u) = C + D < 



1 + 2a - 2a 2 



log 



1 - 


2au \ 


Vi 


+ 4a 


1 - 


2au 




+ 4aJ 



u 2y/l + 4a(2 - a 



(A.2) 



Constants of integration Cq and D should be determined to be a regular function at the 
horizon. So we here choose D = and get 



F (u) =C = C, (const.). 
By using this solution, one can get an equation for F\{u) from ( 14. 61) . 

= (u 2 (l - u)(l + u - au 2 )F'^j + i^~^u(2 + 3u - Aau 2 ^. 



(A.3) 



(A.4) 
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A general solution is 



. - , I + 2a -2a 2 
F 1 (u) =C 1 + D 1 [ - 7^==^ 7 log 



/ l-2a«\ 



2-s/l + 4a(2 - a) 
1 1 



1 — 2au 



. 1 + . 



« 2(2 -a 
— - (iCb - D^j log (l - u) . 



log 1 + u — aw' 



Again, removing the singularity at the horizon, the constant Di should be 

D l = iCb. 

We also impose a boundary condition F±(u = 1) = 0, so as to fix the constant C\, 



Cl -!. 1 . + 2a - 2a ',- : 



. 2 VI + 4a(2 - a) 
Therefore the final form is 



( x l~2a ^ 
VI +4a 
1 -2a 

\ + ym^j 



i + 



2(2 - a) 



log (2 - a 



/••;(,/) = iC6<j - ^ + 20 - | log 



2Vl + 4a(2 - a) 



+ 1 - 



■log 



/ l-2au\ 

VI + 4a 
1 -2a 

V ~ y/T+4a/ 
1 + tt — au 2 



log 



/ 1 - 2au \ 

1 + ; X 

VI + 4a 

1 - 2a 

V + v / iT4^y 



u 2(2 - a) 
A differential equation for Gi(u) is 

= (V(l -u)(l + u-au 2 )G[) 

A general solution is 

5i 



C6 2 « 1 



3a 



2(1 + a) 



— «). 

a / 



(l + 2a - 2a 2 ) (Co 2 + 2£> 1 (1 + a)) 

+ - 7 V Mog 

4 v / TT4^(l + a)(2-a) 



/ 1 - 2aw \ 

yi + 4a 
1 — 2aw 

V 1 + 7r+la7 



Co 2 + 2^(1 + a) / / n 1 / 2 \\ 
2(l + a)(2-a) (^g (l -«j - -log (l+ti-au jj, 
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and the constant Di might be fixed as 

Cb 2 



D 1 



2(1 + a)' 

From the condition G\(u = 1) = 0, we can fix the constant C\ as 

c Cb2 

1 2(1 + a)' 

So we obtain the final form, 



Cb 2 / 1 - 

Gi(u) = — ■ - ( + (A.9) 



2(1 + a) V u 



Appendix B. Perturbative solutions for $ + 

From the equation (14. 10[) . we have a differential equation for Fq(u), 

=((!-«) (!+„-«") (l-j^)^)'. (B.l) 



A general solution is given by 



~ D f 18a(2 -a) 1 - 10a - 2a 2 

F {u) -C - 2 (2-a)3\(i + 4a)(2 + 2a-3au) ~ (l + 4a) 3 / 2 g 



/ l-2a«\ 

VI + 4a 
1 — 2au 

1 + 



+21og(l-u) -log(l+M-au 2 ) I. (B.2) 



Since the function F should be regular at the horizon, we choose D = and get 

F (u) = C = C, (const.). (B.3) 
Substituting the solution to the equation (I4.10p . we get an equation for Fi(u), 
= ( C 1 - u)(l + u - au 2 ) (l - 3Q , u) 2 F[ 



2(1 + a) 



A general solution is given as 



F lM = c 1 + 2 ' 1 + °lV{- 18a < 2 - a 



(2 - a) 3 [ (1 + 4a) (2 + 2a - 3aw) 
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(1 - 10a - 2a 2 ) 
+ (l+4a) 3 / 2 g 

+ log ^1 + u — an 2 j | 

iC(2 -afb- 4L> 1 (1 + a) 2 / \ 
f — — ^r 1 log (1 - uj . 



( \-2au\ 

I 1 ; \ 

yj +4a 
1 — 2aw 



(2 - o)3 

The constant of integration Di should be 



(B.5) 



1_ 6 4(1 + a) 2 ' 

so that the singularity at the horizon would be removed. In addition, we require the 
condition F\{u = 1) = to fix the constant C\, 



d = i- 



Cb f 9a 



1 - 10a - 2a 2 
2 - a \ 1 + 4a 2(1 + 4a) 3 / 2 

Then we get the final form of the solution 

~ . , . Cb f 27a 2 ( l-u 
F i( u ) = K 



log 



( x l~2a ^ 
' " VI + 4a 



1 -2a 



\ 1+ v / TT4a7 



- log | 2 a 



2 - a I 1 + 4a V 2 + 2a - 3au 



1 - 10a - 2a 2 . . 
+ -77^^^77^( log 



2(1 + 4a) 3 / 2 
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+ o log o 
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/ 1 - . \ 

yj + 4a 

i 1 -2a 



-log 



/ l-2au\ 
I 1 + ; \ 

yj + 4a 

i 1 -2a 



(B.6) 



Similarly we have a differential equation for Gi(u), 



= ^1 +u-au 2 ^j{] 



Cb 2 



2(1 + a) 
3a 2 



it v 2(1 + a) /V 2(1 + 
A general solution of this equation can be obtained, 



(B.7) 



Gi(«)=Ci 



9a (#i + 4C(1 + a) 2 6 2 log(3a) 
(2-a) 2 (l +4a) 



2 + 2a — 3aw 
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x log 



/ l-2au\ 

VI +4a 

1 — 2aw 
1 + ; 



C6 2 f(2-a)(14 + 31a + 8a 2 ) + 24(1 + a) 3 log(3a)) + 6£>i(l + a) . N 

^ ; — r '- log (1 -u) 

6(2-a) 3 (l + a) & V / 



4C(l + a) 2 6 2 

C(2-a)(14-21a-84a 2 -76a 3 )6 2 + 6Di(l + a)(l + 4a), / 2 \ 

- — " 12(2 - B )»(i + a)( i + 4.) — " ~ log I 1 + " " au ) 
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' (2 - a) 2 (l + 4a) V2 + 2a - 3a U/ 
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' (2 -a) 3 Ll2(M) 

2C(1 + a) 2 o 2 
"(2-a) 3 (l + 4a) 3 / 2 

xi (l-10a-2a 2 + 



log -a log ^1 — -a) 



(1 + 4a) 3 / 2 ) 



x(log(3a M )log^l- r _^__ ; 
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2au \ 



+ Li 2 



2au \ 
1 - VYTIa) 



x ( log (Sau) .og (l - YT ^ m ) + Li 2 ( TT ^L=) ) } ■ 



(B.8) 



The constant of integration Di might be fixed to remove the singularity u — 1, 

(7o 2 



#1 = 



(2 - o)(14 + 31a + 8a 2 ) + 24(1 + a) 3 log(3a) . 



6(1 + a) 

Another constant of integration C 1 is used to satisfy the condition Gi(u = 1) = 0. The 
final expression of the solution is 



G 1 (u) = Cb 2 



9a 2 (14 + 31a + 8a 2 ) / 1 - u 
'2(l + a)(l + 4a)(2-a) 2 V2 + 2a-3aw, 

(1 + a) (3a(2 - a) (5 + 2a) - 2(1 + a)(l - 10a - 2a 2 ) log(3a)) 

(2-a) 3 (l + 4a) 3 / 2 
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\ 1 - VT+~W 



+Li 2 



2aw 



Li 5 



2a 



.1 - \/T+4a/ V 1 - VI + 4a 

(l-10a-2a 2 -(l + 4a) 3/2 



x ( log u log ( 1 



2au 



1 + VI + 4a 



+Li 2 



2aw 



1 + VI + 4a 



— Lio 



+ log(3a) log 



2a 



/ 2au \ 

1 + VI + 4a 



2a 



1 + VI + 4a 



\ 1 + vT+W 

(B.9) 



Appendix C. Perturbative solutions for h xy 

From the equation of motion (15. 2p . one can get an equation for F (u) 

(-(1 - u)(l + u - au 2 )F^j' = 0. 
A general solution is given by 



(C.l) 



F {u) = C + 



Dn 



2- a\ 2 VI + 4a 



log 



/ l-2aw\ 

VI + 4a 
1 — 2aw 
\ 1 + ; / 

\ VTTW 
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-log(l-u) + ^log(l + «-a« 2 )J. (C.2) 



Constants of integration Co and D should be determined to be a regular function at the 
horizon. So we here choose D = and get 



F (u) =C = C, (const.). 

By using this solution, one can get an equation for Fi(u), 

f-(l - u)(l + u - au 2 )F[) = ( 1 + au 2 ). 

\u J 2 — a u 2 V / 

A general solution is 

Fi(u) =d- ]_ 2 (tC6 + (2 - o)£>!) 



(C.3) 



(C.4) 



x 



log 



/ l-2au\ 
VI + 4a 



1 + 



1 — 2aw 



^—t(iCb(\ - a) - (2 - a)Di) log (l - uj. 



log ^1 + -u — a-u 2 j | 



(2 - a)' 

Removing the singularity at the horizon, the constant Di should be 

D 1 = iC l —b. 
2- a 

We also impose a boundary condition Fi(u = 1) = 0, so as to fix the constant Ci, 

Cb 



(C.5) 



( x l-2a \ 



1 ~ 2(2 -a) \7r+4^ 
Therefore the final form is 



log 



VI +4a 
l I -2a 



log (2-a)|. 



Fx(u) = i- 



C6 



2(2 - a) \ VTT4^ 



log 



/ 1 - 2au \ 
l-2a 



log 



+ log 



V VTTw 

1 + u — aw 2 



/ 1 - 2a« \ 
/ 1 + ; 1 

VTT4^ 

i 1 -2a 



2 - a 

A differential equation for Gi(-u) is 



(^(l-uXl + u-au 2 )^)' 
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Cb 2 



(C.6) 



(C.7) 



The equation gives us the solution 

C(l - 2a)b 2 + 3Dt 



GAu) = d 



2Vl + 4a(2 - a) 



log 



1 - 


2 aw 


Vi 


+ 4a 


1 - 


2 aw 



\ vT+4a/ 



Co 2 - Di 



2(2 - a 

and the constant Di might be fixed as 



2 log ( 1 — u ) — log (l + ii — <xu 



D 1 = Cb 2 . 



We can also fix the constant C\ as 



c i = log 



A l-2a \ 



VI + 4a 



1 - 2a 



. 1+ , . 



Then we obtain the result, 

Cb 2 



GAu) 



VI + 4a 



log 



/ l-2au\ 

VI + 4a 
1 -2a 

V ~ VTT4^/ 



log 



/, l-2aw\ 
/ 1 + , \ 

VI +4a 



1 - 2a 



. 1+ , , 



(C.8) 



(C.9) 
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